Quantum critical behavior and trap-size scaling of trapped bosons in a 
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We study the quantum (zero-temperature) critical behaviors of confined particle systems described 
by the one-dimensional (ID) Bose-Hubbard model in the presence of a confining potential, at the 
Mott insulator to superfluid transitions, and within the gapless superfluid phase. Specifically, we 
consider the hard-core limit of the model, which allows us to study the effects of the confining 
potential by exact and very accurate numerical results. We analyze the quantum critical behaviors in 
the large trap-size limit within the framework of the trap-size scaling (TSS) theory, which introduces 
a new trap exponent 9 to describe the dependence on the trap size. This study is relevant for 
experiments of confined quasi ID cold atom systems in optical lattices. 

At the low-density Mott transition TSS can be shown analytically within the spinless fermion 
representation of the hard-core limit. The trap-size dependence turns out to be more subtle in 
the other critical regions, when the corresponding homogeneous system has a nonzero filling /, 
showing an infinite number of level crossings of the lowest states when increasing the trap size. At 
the n = 1 Mott transition this gives rise to a modulated TSS: the TSS is still controlled by the 
trap-size exponent 6, but it gets modulated by periodic functions of the trap size. Modulations of 
the asymptotic power-law behavior is also found in the gapless superfluid region, with additional 
multiscaling behaviors. 
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I. INTRODUCTION 

The impressive progress in the experimental manipula- 
tion of cold atoms in optical lattices (see, e.g., Ref. [1] and 
references therein) have provided a great opportunity to 
investigate the interplay between quantum and statisti- 
cal behaviors in particle systems. Cold atoms in optical 
lattices can be used to study many-body phenomena in 
dilute gases, such as quantum Mott- Hubbard transitions 
for bosonic atoms, see, e.g., Refs. [2-9]. An important 
feature of these experiments is the presence of a confining 
potential which traps the particles within a limited spa- 
tial region of the optical lattice created by laser-induced 
standing waves. The theoretical framework [10] is based 
on the Bose-Hubbard (BH) model [11] in the presence of 
a confining potential coupled to the particle density, i.e., 
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where (ij) is the set of nearest-neighbor sites, Ui 
is the particle density operator. 

We consider a power-law trapping potential 



V{r) 



(1) 



(2) 



the particles are trapped. The inhomogeneity due to the 
trapping potential strongly affects the phenomenology of 
quantum transitions in homogeneous systems. 

The homogeneous BH model without trap undergoes 
Mott insulator to superfluid quantum transitions driven 
by the chemical potential /z, whose low-energy proper- 
ties are described by a nonrelativistic U(l)-symmetric 
bosonic field theory [11], which is characterized by the 
dynamic exponent z = 2.^ In the presence of a confin- 
ing potential, theoretical and experimental results have 
shown the coexistence of Mott insulator and superfluid 
regions when varying the total occupancy of the lattice, 
see, e.g., Refs. [7, 10, 13-22]. However, at fixed trap size, 
the system does not develop a critical behavior with a 
diverging length scale [13, 17]. 

Criticality can be recovered only in the limit of large 
trap size. As put forward in Refs. [23, 24], this critical 
regime can be described in the framework of the trap- 
size scaling (TSS) theory, where the critical behavior is 
cast in the form a TSS with a nontrivial trap exponent 
6, which determines how the length scale of the critical 
modes at the critical point diverges with increasing trap 
size, i.e., ^ ~ Z^. For example, let us consider a standard 
scenario (see, e.g., Ref. [25]), in which the quantum T — 
transition of the the homogeneous d-dimensional system 
has one relevant parameter fi, with critical value fic- The 



where r = \x\, v and p are positive constants and I = l/v 
is the trap size. Experiments are usually set up with 
a harmonic potential, i.e., p — 2. Far from the origin 
the potential V{r) diverges, therefore (n^) vanishes and 



The special transitions at fixed integer density belong to a dif- 
ferent universality class [11], described by a relativistic U(l)- 
symmetric bosonic field theory, which is the {d + l)-dimensional 
XY universality class [12]. Thus its dynamic exponent is z = 1. 
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simplest TSS Ansatz [23] for the asymptotic behavior 
of the free-energy density in the presence of a confining 
potential (2) is 

F{fi, T, I, x) = r^('^+") xi-^), (3) 

where x is the distance from the middle of the trap, j2 = 
H — He z is the dynamic exponent and v = l/Un where 
is the renormalization-group (RG) dimension of /i. 
Moreover, any low-energy scale at T = 0, and specifically 
the gap, is expected to behave as 

A = l-^^'Viiif/"). (4) 

The above TSS has been verified by analytical and ac- 
curate numerical calculations [23] within the quantum 
XY chain in a space-dependent transverse external field, 
which acts as a trap for the spinless fermions of its 
quadratic Hamiltonian representation, which can be ob- 
tained by a Jordan- Wigner transformation. 

The general features of the TSS at the Mott to super- 
fluid transitions in rf-dimcnsional BH models have been 
discussed in Ref. [23]. Beside the dynamic critical ex- 
ponent z and the RG dimension of fjb, which control 
the critical behavior of the homogeneous system and can 
be determined from the corresponding continuum the- 
ory [11, 25], the TSS requires the trap exponent ^, which 
can be derived by a RG analysis of the corresponding per- 
turbation [23]. For one- and two-dimensional systems, 
the critical exponents entering the scaling formulas (3) 
and (4) are 

z = 2, y^ = l/y = 2, e=p/{p + 2), (5) 

where p is the power of the confining potential (2). 

In this paper we investigate the quantum critical be- 
haviors of the one-dimensional (ID) BH model at zero 
temperature in the presence of a confining potential, at 
the Mott to superfluid transitions, and within the gapless 
superfluid phase. The ID BH model in the presence of a 
confining potential is of experimental relevance in optical 
lattices, where quasi ID confined particle systems have 
been realized, see, e.g., Refs. [1, 4, 5, 9]. 

Specifically, we consider the hard-core limit, f/ ^ oo, 
of the ID BH model, which allows us to study the ef- 
fects of the confining potential by exact and very accu- 
rate numerical results. The hard-core limit implies that 
the particle number rii per site is restricted to the values 
Ui = 0,1. In this limit the model can be mapped into 
the XX chain model with a space-dependent transverse 
external field, 

i 

- Y.[i, + V{x,)]Sl (6) 

i 

where Sf = (Tf/2 and cj" are the Pauli matrices, which 
are related to the boson operators bi by = + bi, 



af = i{bl — bi), erf = 1 — 2&]6i. In the following we fix 
J = 1. Then, by a Jordan- Wigner transformation, one 
can further map it into a model of spinless fermions, see, 
e.g., Ref. [25]. 

In the absence of the trap, the ID hard-core BH model 
has three phases: two Mott insulator phases, for /z > 1 
with {rii) ~ and for /i < — 1 with (n^) = 1, separated 
by a gapless superfluid phase for ]/i] < 1. Therefore, 
there are two Mott insulator to superfluid transitions at 
/Lt = 1 and /Lt = — 1. At both transitions the exponents 
controlling the critical behavior are those reported in Eq. 
(5) . The gapless superfluid phase is instead described by 
a free massless bosonic field theory with dynamic expo- 
nent z = 1, see, e.g., Ref. [26]. 

The effects of the confining potential at the low-density 
transition {fj, = 1) has been already studied in Ref. [23], 
where TSS has been shown to emerge by analytical cal- 
culations exploiting the spinless fermion formulation of 
the ID hard-core BH model. 

In this paper we extend the study of the trap-size de- 
pendence to all critical regions of its phase diagram. We 
present results for the trap-size dependence in the gap- 
less superfluid region and at the n = 1 Mott insulator to 
superfluid transition. 

For this purpose we exploit the free spinless fermion 
representation of the ID hard-core BH model, which al- 
lows us to perform computations for very large systems, 
since they only require the diagonalization of a i x i 
matrix where L is the number of lattice sites. We obtain 
numerical results for chains of size L, with a trap of size 
/ centered in the middle of the chain (unless explicitly 
stated, we consider odd Ls so that the middle of the trap 
coincides with the middle site of the chain); we choose 
L large enough to have negligible finite-L effects; we are 
able to obtain results correct to machine precision for I 
up to 0(10^ ). We then analyze the quantum critical be- 
haviors in the presence of the trap within the framework 
of the TSS theory. 

The trap-size dependence shows subtle effects in the 
parameter region where the homogeneous model without 
trap has a nonzero filling /, i.e., for /z < 1, therefore in 
the superfluid region and at the n = 1 Mott transition. 
This is essentially related to the presence of level cross- 
ings at finite trap size. They arise because the particle 
number is conserved, i.e., the particle number operator 
N = rii commutes with the BH Hamiltonian (1) even 
in the presence of the trapping potential; thus the eigen- 
vectors do not depend on /i, even though the eigenvalues 
do. In the presence of the trapping potential (2), the 
particle number N = {N) is finite and increases as N I 
with increasing the trap size /. Therefore, as / — > oo, 
there is an infinite number of ground-state level cross- 
ings where A'' jumps by 1 and the gap A vanishes. 

As we shall see, this phenomenon gives rise to a new in- 
teresting scenario at the n = 1 Mott transition, requiring 
a revision of the TSS Ansatz (3) and (4) into a modulated 
TSS: the TSS is still controlled by the trap-size exponent 
= p/{p + 2), as in the case of the low-density Mott 
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transition, but it gets modulated by periodic functions 
of the trap size. 

We provide numerical evidence of universality at the 

low-density and n = 1 Mott transitions by considering 
a more general ID hard-core BH model with nearest- 
neighbor density-density interactions, corresponding to 
the so-called XXZ model. Since the corresponding 
fermion representation is not longer quadratic, the nu- 
merical results are obtained by using density matrix 
renormalization group (DMRG) methods. 

We also show that the trap-size dependence in the 
gapless superfluid phase is characterized by power-law 
asymptotic behaviors which are modulated by periodic 
functions of the trap size. Moreover, it shows a multiscale 
behavior characterized by different length scales diverg- 
ing with different power laws with increasing trap size, 
associated with the smooth modes and with the modes 
at the Fermi momentum kp ~ irf. 

The paper is organized as follows. In Sec. II we dis- 
cuss the TSS at the low-density Mott transition in the ID 
hard-core BH model. We extend the analysis of Ref. [23] , 
providing the TSS of several observables and checking its 
universality within the XXZ model. In Sec. Ill we discuss 
the local density approximation to determine the parti- 
cle density in the presence of the confining potential. In 
Sec. IV we study the trap-size dependence within the 
gapless superfluid phase, presenting results for the gap, 
the particle density and its correlators, the one-particle 
density matrix and the von Neumann entanglement en- 
tropy. In Sec. V we consider the n = 1 Mott insulator 
to superfluid transition at fj, = —1. We present results 
at fixed trap size for the XX and XXZ models, and show 
that their trap-size dependence is described by a modu- 
lated TSS. Finally, in Sec. VI we draw our conclusions. 
In App. A we report some details on the numerical cal- 
culations presented in paper. App. B reports some re- 
sults for the homogeneous ID hard-core BH model with 
open boimdary conditions, showing that modulated scal- 
ing behaviors of the gap, the particle density, and the 
subleading corrections of the entanglement entropy, are 
already present in the finite size behavior of the homo- 
geneous model within the superfluid region; we provide 
exact formulae for the total particle number and for the 
particle density at the middle of the chain. 



II. TRAP-SIZE SCALING AT THE 
LOW-DENSITY MOTT TRANSITION 



A. The TSS limit 

In order to show the existence of a nontrivial TSS limit 
around fx = 1, i.e., at the transition between a low- 
density superfluid and the empty vacuum state (which 
may be named n = Mott phase), wc exploit the ex- 
act mapping of the XX model into a model of spinless 



fermions, by the Jordan- Wigner transformation 



(7) 



j<i 



ay = il[{l-2c]cj){cl-c,), 

j<i 

af = 1- 2c] c. 



which leads to the Hamiltonian 
He = c\hijCj, 



hij = Sij - ^Sij-i - ^Sij+i + + V{xi)]6i. 



(8) 



with /i = /i — 1. In the fermion representation the Hamil- 
tonian can be easily diagonalized by introducing new 
canonical fermionic variables r]k = <pkiCi, where ^ sat- 
isfles the equation 



obtaining 



hij(j)kj = (^k<i>ki, 



He = y^^0Jkvlr]k- 



(9) 



(10) 



The ground state contains all Ty-fermions with ojk < 0; 
the number of filled energy levels is N. The energy gap 
is 



A = min \uik\ 

k 



(11) 



Since (j)ki is an orthogonal matrix, the expectation values 
of the c-operators can be obtained by using the inverse 
relation Cj = J2k 'PkiVk- 

The above equations have a nontrivial TSS limit 
around 



/X = /It — 1 = 0, 



(12) 



i.e., at the transition between a low-density superfluid 
and the empty vacuum state. We consider the contin- 
uum limit of Eq. (9), by rewriting the discrete differ- 
ences in terms of derivative expansions. Near the critical 
point p, = and for sufficiently small values of k (this 
is required by the smoothness hypothesis underlying the 
continuum limit), we obtain 



M+(W- 



1 rf2 



1 



+ ... 



2 24 dx'^ 

where (j)k{x) = <pkx- Then, by replacing 



(f)k{x) = uJk(l>k{x), 
(13) 



(14) 
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and neglecting terms which are suppressed in the large-Z 
hmit, we obtain 

(15) 

where (pk{X) = note that flk is indepen- 

dent of fir- Recalling that z ~ 2 and = 2 at this tran- 
sition, we infer 6 ^ p/{2+ p), as also obtained by RG 
arguments [23]. Therefore, the solutions of Eq. (15) de- 
termine the TSS in the limit /x = /Lt — 1 — > and I — >■ cx), 
keeping fir = l^^p- and X = l~^x fixed. On the other 
hand, subleading terms, such as the d'^/dx'^ term in the 
l.h.s. of Eq. (13), give rise to 0{l~^^) corrections in the 
TSS hmit (14) of the Eq. (9). Thus, we expect that the 
TSS of any quantity is approached with 0{l^^^) scaling 
corrections. 

For p = 2, we obtain 

l]fe = rjfe - /ir = 2^/^(fc + 1/2), fc>0, (16) 

^^{X) = exp(-XVV2), 



7ri/42'=/2(fc!)i/2 
where X = x/l^^'^ and H^ix) are Hermite's polynomials 



(-1)* 



/dx^ 



Hk{x) 

For p = 4, we can solve numerically Eq. (15) by Nu- 
merov's method; the resulting energy levels are 0,q — 
0.667986, ni = 2.39364, ^2 = 4.69680, fl^ = 7.33573, 
^4 = 10.2443, ^5 = 13.3793, etc.; Bohr-Sommerfield 
quantization formula gives 

« [74(fc + l/2)f 
74 = V^r(7/4)/r(5/4) 9^ 1.27082, (17) 

which is accurate to 0.1% already for Ct^. For large k and 
generic p, the semiclassical limit gives 



(18) 



For p — )■ oo, Eq. (15) becomes equivalent to the 
Schrodinger equation of a free particle in a box of size 
L = 21 with boundary conditions ipi—i) = <p(l) — 0, 
leading to 



fc > 0, 



(19) 



MX) = sin[-ik + l)iX + l) 



where X = x/l. 
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FIG. 1: The rescaled energy gap / A (below) and the rescaled 
particle density {no) in the middle of the trap (above) for 
p — 2 {6 = 1/2), compared with the predictions of TSS. We 
report results for several values of p, and / > 10. They are 
perfectly consistent with the analytical TSS. 
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B. TSS of observables 




1. The energy gap 

The existence of the TSS limit implies that any low en- 
ergy scale behaves as l"'^^ £ (fir) , where £{nr) is a scaling 



FIG. 2: The rescaled energy gap Z^" A (below) and the rescaled 
particle density {no) in the middle of the trap (above) for 
p — 4 (9 = 3/4), compared with the predictions of TSS. We 
report results for several values of jl and I > 10. 
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FIG. 3: The rescaled energy gap l^^ A (below) and the rescaled 
particle density (no) in the middle of the trap (above) for 
p — >■ oo (S = 1), compared with the predictions of TSS. We re- 
port results for several values of Jl and I > 10. They approach 
the analytical TSS results with increasing /. 



function. Specifically, in the case of the gap IS. = Ei- 
we have 



- En 



k 



(20) 



which can be easily computed from the solution of Eq. 
(15). The results are shown in Figs. 1,2, and 3, respec- 
tively for p = 2, 4 and in the limit p — oo. For any p, 
the scaling function SAifJ-r) shows a triangular structure 
for /ir < and it is linear for /i^ > — f^o < 0- In Figs. 
1, 2, and 3 we also show results obtained by numerical 
diagonalization at fixed /. They clearly approach the an- 
alytical TSS computations. Corrections to scaling turn 
out to be very small for p — 2 and p — A. 



This is confirmed by the analytical results of the previous 
section, which lead to 



k: nk<o 



(23) 



where we used the fact that {rf^rik) = 1 if J7fc < and 
zero otherwise. ipk{X) are the normalized eigenfunctions 
of Eq. (15). Note that T) depends on /i^ only through 
the number of negative energy levels N , i.e., the number 
of levels included in the sum of Eq. (23); therefore, it 
vanishes for Hr > —^o and it falls on a discrete set of 
curves as a function of X (with jumps at VL^ = 0, i.e., 
zeroes of A). Since '^k{X) = (— l)'^(/?fe(— X), only even 
fcs contribute to 2?(/Xr, 0). 

Specifically, for p = 2 we obtain 



p(,.,o)=2vv-v.|:[(^-^--i)'']^ 



Jmax — max 



m 

(O, [|Ai.|/y8-l/4 



(24) 



where \x\ = floor(x) is the largest integer not greater 
than X. For p cx) we have 



2?(Mr,0) 



0, 



V2H/7r + 1/2 



(25) 



Numerical results for the particle density at the origin 
are shown in Figs. 1, 2, and 3, for p = 2, 4, oo respectively; 
they fully support TSS. Note the peculiar plateaus and 
the discontinuities in the particle density at negative val- 
ues of the scaling variable /i^- Moreover, asymptotically 
for /ir — > — oo, {tiq) ^ ImT^^, which matches the criti- 
cal behavior for /2 < in the absence of the trap [25]. 
Numerical results for {ux), showing the dependence on 
the distance form the middle of the trap, are shown in 
Fig. 4, for -0.01 < < and ? > 10. They show that 
the quantity {ux) approaches the analytical functions 
obtained using Eq. (23). 



3. The particle density correlator 



2. The particle density 

We now consider the expectation value and correlators 
of the particle density 

Tlx = blbx = = clcx = ^ (j}*k, x'l^k2xVtVk2- 

kl ,k2 

(21) 

(ux) is zero for /i > 0, but it can be nonzero for fi < 0. 
Since the RG dimension of the particle density is given 
by Un = d + z - = 1, we expect 

(rix) ^l-''V{iir,X). (22) 



Let us consider the particle density correlation 

Gn{x) = {noTix) - {no){nx) = ^(cto^^)c- (26) 

Like the particle density, it vanishes for p. > 0, and it can 
be nonzero for /2 < 0. In the TSS limit we expect the 
scaling behavior 

Gn{x)^l-^'gn{flr,X). (27) 

Straightforward calculations, i.e., writing G„ in terms of 
77-operators and then using the Wick theorem to compute 
the resulting ground-state expectation values, show that 
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FIG. 4: The rescaled density f iux) iox p — 2 (above), p = 4 
(middle), and p — > oo (below), compared with Eq. (23). A'^ 
is the number of negative energy levels with Q.^ < 0. Here 
and in the following, set of parameters {fi, I) are differentiated 
by symbol color (matching A'^), shape, and filling color. The 
results approach the analytical A''-dependent TSS results, cf. 
Eq. (23), with increasing /. 



the scaling function tj„(/ir, X) can be written in terms of 
the eigensolutions of Eq. (15): 

gn{i^r,X)= (28) 




FIG. 5: The rescaled connected correlation l'^^Gn{x) for p = 
2, I > 10, /i > 0.95, and A*' < 6. The results approach 
the analytical A^-dependent TSS results, cf. Eq. (29), with 
increasing /. 



then, using the completeness relation ipk(X)ipk{Q) = 
5{X), we obtain 



Qn{t^r-,X) 



nfc<0 



Vk{X)<Pk{^) 



(29) 



note that odd fcs do not contribute to the sum. Like 
the particle density, Qn{^J'^r,X) is nonzero only for < 
—Clo < 0, and depends on /i^ only through N (actually, 
only through [{N + 1)/2J). 

From the numerical data, we find that P^Gn{x) as a 
function of X approaches rapidly GnifJ-n X). The results 
for p = 2 are shown in Fig. 5. 



4- The one-particle density matrix 

We now consider the one-particle density matrix de- 
fined as 



Gb{xt,Xj) = {bib J 



(30) 



where ct± = {a'= ± i<jy)/2. For A > 0, since the ground 
state is empty, Gb{xi,Xj) = 0. But, like the particle den- 
sity (note that Gb{x,x) = {fix)), it can be nonzero for 
jl < 0. Its scaling behavior is also determined by the RG 
dimension yi, of the field associated with the boson oper- 
ator bi, which is [25] yb = 1/2. Considering specifically 
the correlations with one boson operator at the origin, 
i.e., Gb{x,0), we expect 



Gb{x,0) w l-''M{^lr,X). 



(31) 



Numerical results can be obtained using the method of 
Ref. [27]. They confirm the scaling behavior (31), i.e., 
l^Gb{x,0) appears to approach a function of fj,r and X 
in the large-trap limit. Fig. 6 shows results for p — 2. 
Again, the dependence on fj.^ is only through the number 
N of negative energy levels. 
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FIG. 6: The rescaled correlation l^Gb{x,0) for p = 2, 
I > W, ^ > 0.95, and < 6. With increasing /, they converge 
toward A^-dependent TSS curves. 
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FIG. 7: The half-lattice von Neumann entanglement entropy 
Si/2, cf. Eq. (32), for p = 2. 



5. The von Neumann entanglement entropy 

We consider the von Neumann entanglement entropy 
S{Ia]L) in the presence of the confining potential, de- 
fined by dividing the chain in tw^o parts of length I a and 
L — Ia- The entanglement entropy trivially vanishes for 
/I > 0, but it can be nonzero for /l < [28]. It can 
be computed using the techniques of Ref. [29] . We con- 
sider the half-lattice von Neumann entanglement entropy 
S{L/2]L) for even L and open boundary conditions in 
the presence of the trap of size I (in this case the trap is 
centered between the two central sites of the chains) . Its 
large-L limit, 



5i/2 = hm S{L/2;L) 



(32) 



depends on the trap size I only. The TSS limit of 6*1/2 
turns out to depend only on TV, and therefore it is a 
function of the scaling quantity fir, as shown in Fig. 7 
for the harmonic potential. 



C. Universality of the TSS 

In order to check the universality of the TSS, we con- 
sider the XXZ chain model 

Hxxz = -^{Si S^^ 1 + Sf Sf^ i+jz S- S-^ 1 ) 

i=l 

-J2[f, + Vix.)]S^. (33) 

i=l 

In terms of bosonic operators bi, the jz term corresponds 
to nearest-neighbor density-density interactions, i.e.. 



J 



,,;^(n,-l/2)(n,+i-l/2). 



(34) 



where iii = h\hi. 

For < 1, the homogeneous XXZ chain model (i.e., 
with V = 0) undergoes two Mott insulator to superfluid 
transitions at /i^ = ±(1 -I- j^) in the same universality 
class of those of the XX model at /i = ±1. Indeed, the jz 
term is irrelevant at these transitions: the RG dimension 
of the coupling jz is yj^ = — 1 [25]. We therefore expect 
that the jz term remains irrelevant also in the TSS limit. 
The main difference concerns the dominant scaling cor- 
rections, which are expected to be 0(f"^), thus 0(r''), 
when this term is present. Therefore they are expected 
to be larger than those of the XX model, where they are 
0(/-2«). 

The Hamiltonian (33) is no longer equivalent to a 
quadratic ferniionic Hamiltonian. The properties of the 
model can be studied numerically via DMRG. Results for 
the gap and the particle density in the middle of the trap 
at the low-density Mott transition, around /i^ (we define 
p, = fi — /i^), are shown in Fig. 8 for p = 2. They clearly 
support universality of TSS. As expected, corrections to 
scaling are larger for jz 7^ than for jz — 0. 

It is interesting to notice that energy differences and 
expectation values over states with 2;-component of the 
total spin M ^ L/2 and M = i/2 — 1 depend on /i and 
jz only through This implies that, in the scaling 
region, vacuum expectation values (as functions of fir) 
are independent on jz for /i^ <^ — ^^2; the energy gap A is 
independent on for /i^ > — (Hi -|-n2)/2. This is clearly 
shown in Fig. 8. 



III. THE LOCAL DENSITY APPROXIMATION 
OF THE PARTICLE DENSITY 

The homogeneous ID hard-core BH model, i.e., the 
model (1) for J7 — > 00 and V{x) = 0, has a nonzero filling 
below the low-density Mott transition, i.e., for fi < 1. In 
the infinite-chain limit L — ?► 00, the filling / for < 1 
is given by [25] 



N 



1 



■ arccos /i. 



fi — cosnf. (35) 
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FIG. 9: The particle density vs. x/l for p = 2 and several 
values of /i. We show results from the LDA, cf. Eq. (37), and 
from numerical calculations on a large chain with I = 200. 
The differences are hardly visible in the figure. 



when ficS < — 1, for 

x/l<i-l-fiy/P, (38) 
and a vanishing particle density when > 1, for 

x/l> {1- fi^/P. (39) 



FIG. 8: Scaling plot of the energy gap (below) and of the 
particle density at the origin (above) for the XXZ model with 
p = 2 and = -1/4, 1/2. The data at = -1/4 and 1/2 
approach the TSS functions, with increasing from opposite 
sides. 

The corresponding Fermi momentum is kp — nf. For 
^ < — 1, / = 1 independently of ^. In the following, / 
will always denote the value for the infinite homogeneous 
chain. 

In the presence of a space-dependent confining poten- 
tial, the so-called local density approximation (LDA) es- 
timates the spatial dependence of the particle density by 
taking the value of the particle density of the homoge- 
neous system at the effective chemical potential 

fi,s{x) = fi + {x/l)P. (36) 

The LDA has been widely used to get quantitative infor- 
mation on the behavior of BH models in a confining po- 
tential, and, more generally, of inhomogeneous systems, 
see, e.g., Refs. [13, 16, 17, 30, 31]. 

The LDA of the particle density reads 

("-x)lda = /3ida(a;/0 = (37) 

for fics{x) > 1, 
(I/tt) arccos/ioff(a;) for — 1 < Atefi(a;) < 1, 

1 for Hes{x) < -1. 



In Fig. 9 we compare the LDA of the particle density with 
numerical results for p = 2 and the trap size I = 200. 
Note the flat regions related to the n = 0, 1 Mott phases, 
already observed in experimental and numerical works, 
see, e.g., Refs. [7, 10, 13]. Analogous results are found 
for other powers of the confining potential. The LDA 
provides a good approximation of the particle density, 
which improves with increasing trap size. The differences 
of the trap-size dependence from the LDA results show 
a nontrivial scaling behavior; they will be considered in 
the next sections. 

Using Eq. (37) we can also obtain the LDA of the total 
particle number: 

POO 

Mda = V'("-x)ida ^21 pidaiy) dy = C{fi)l, (40) 

where c(/i) is a finite function of fi, which can be easily 
computed by integrating pida(?/)- Comparing with the 
numerical results, we find excellent agreement, i.e., 

\c{fi)l - N,^,,t\ < 1 (41) 

(note that the A'oxact is an integer number, while the 
LDA is a continuous linear function). Eq. (41) implies 
that, in the large-^ limit at fixed fi, the total particle 
number increases as 



This would imply the presence of a plateau at n = 1 



(42) 
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IV. TRAP-SIZE DEPENDENCE IN THE 
SUPERFLUID PHASE 

Wc now discuss the trap-size dependence in the gapless 
superfluid phase for < 1. In the continuum hmit, the 
gapless superfluid phase of the homogeneous system is 
described by a free massless bosonic field theory with 
dynamic exponent z = 1, corresponding to a conformal 
field theory with central charge c = 1, see, e.g., Ref. [26]. 



A. Level crossings in the presence of the trap 

In the gapless superfluid phase, and more generally 
for fi = n — 1 < 0, the ground state contains all the ry- 
fermions with cok < 0, cf. Eq. (10). In the presence of the 
trapping potential (2), level crossings of the lowest states 
occur in the jl-l plane separating the regions with N = k 
and N = k + 1. This is essentially related to the fact 
that the particle number is conserved, i.e., the particle 
number operator N = ^ ■ blbi commutes with the BH 
Hamiltonian (1); thus the eigenvectors do not depend on 
/Lt, even though the eigenvalues do. Since, for < 1 and 
in the absence of the trap potential, the ground state 
has a finite density N/L > 0, and in the presence of the 
trap N is finite and increases as TV ~ the lowest states 
show an infinite number of level crossings as ? — >■ oo (after 
L oo) where the gap A vanishes. Note that the hard- 
core limit, [/ — >■ oo in Eq. (1), does not play any special 
role, so we expect that level crossings at finite trap size 
are a general feature of the BH model in the presence 
of a confining potential, when the homogeneous limit of 
infinite trap size has a finite particle density. 

In the following of this section we present results for 
the trap-size dependence of several observables in the su- 
perfluid phase, for < 1. As wc shall sec, the above- 
mentioned level crossings and the competition of smooth 
modes and modes at the Fermi momentum kp = T^f give 
rise to peculiar modulated trap-size dependencies. 



B. Analytical results for the energy gap at small 

We can infer some information on the trap-size de- 
pendence of the energy gap at /i < and small by 
using the analytical calculations of Sec. II A. Let us in- 
troduce a few definitions. Eq'^ (l) is the lowest energy in 

(k) 

the sector with N = k. Iq is the value of / such that 
E^q''\i) = E^''~^^\l); Eq. (10) implies that this is also the 

doubly-degenerate ground-state energy. The gap shows 
peaks at with ^'^^ < < 4'^'^ ^tik is the value 
of I such that eI^''~^\i) — eIj''^^\i); this is the energy 
of the doubly-degenerate first excited state in the sector 

(k) 

with N = k, where the energy of the ground state is Eq . 
For p = 2 Eq. (16) tells us that, for large trap size, the 



gap vanishes at 



'0 



V2 



(fc + l/2). 



(43) 



and therefore the interval between two zeroes is constant, 
i.e., 



,(fe+i) 



(fe) 



V2 



Moreover, the gap A = minj, \oJk\ has peaks at 



;(fe) 
peak 



where the gap decreases as 



r 



(44) 



(45) 



(46) 



The product between the trap size and the gap, i.e., ^A, 
has a periodic asymptotic behavior. Indeed, defining 



I -I 







(fc+i) /(fc) 



for 1^0^ <l<l, 



(fe+i) 



(47) 



"O "0 

thus < < 1, we obtain a simple triangle- like form 

lA = at{(l)), (48) 
where t{^) is the triangle function 

t{<f>) = 1/2 -\cp- 1/2| (49) 

and a = \/2. Note that l^J^,]^ corresponds asymptotically 
to the value (p — 1/2. 

More generally, for any p, the asymptotic behavior (18) 
for large k of the eigensolutions of Eq. (15) implies 

p + 2 



lik) 
'0 



m 



p = 



2p 



(50) 



Thus, for any power p. the location of the zeroes increases 
linearly with fc, and therefore the intervals between sub- 
sequent zeroes i'q'^^'' — Iq'^ approaches a constant as in 
the p = 2 case. This implies again that the peak of the 
gap decreases as and that the product lA is a peri- 
odic function of I, and, specifically, of the corresponding 
phase-like variable <j), cf. Eq. (47). Moreover, asymp- 
totically, lA must have a triangle-like form for any p, 
analogous to Eq. (48). Notable values are 



(fe) 



,{k) 
'o 



«74|M|-3/*(fc + l/2) (p = 4), 
7r2-3/2|/i|-i/2(fc + i) (p = oo), (51) 



and 



11/4 



7r./W2 



{P = 4), 
(p = oo). 



(52) 



The above results suggest that the gap vanishes with 
a global power scaling A ~ However, its amplitude 
in not a constant, but a periodic function of the trap 
size. This scenario will be confirmed by the results from 
numerical diagonalization at fixed ji < 1. 
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TABLE I: Numerical results in the superfluid region for the 
XX model and several values of n and p, for the asymptotic 
interval Q* between level crossings, the constant a entering 
the asymptotic behavior (56) of the gap, the constant b ap- 
pearing in Eq. (59) of the particle density in the middle of 
trap, and the amplitude of the entanglement length scale Ue, 
cf. Eq. (70). The reported estimates are obtained by extrap- 
olating data at finite trap size; their uncertainty is at most 
one on the last figure. 





P 


Q! 


a 


b 


Be 


-1/2 


2 


0.8265187 


1.03010681 


0.378618 


2.64118 


-1/2 


4 


0.8021659 


1.19898955 


0.440691 


2.26916 





2 


1.3110287 


1.198140234 


0.3813798 


2.622058 





4 


1.1635926 


1.349953898 


0.4297036 


2.327185 


1/2 


2 


2.7337528 


1.3177760 


0.484352 


2.06462 


1/2 


4 


2.0551051 


1.29539790 


0.476126 


2.10028 



C. Results from numerical diagonalization 

We now present numerical results at fixed I and fx, 
obtained by numerical diagonalization of the quadratic 
Hamiltonian (8), see App. A. The results at fixed trap 
size are essentially correct up to machine precision. 



1. Interval between level crossings 

The numerical results show that the interval between 
two level crossing approaches a constant value in the large 
trap-size limit. They turn out to fit the simple Ansatz 



(fe+i) 



lik) 



Qi +C2I' 



(53) 



for sufficiently large trap size, ^ ^ 10, which allows us to 
get accurate estimates of the large-Z limit Q'^ . We report 
a selection of results in Table I. 

Note that the results for p = 2 are quite close to, and 
clearly approach, the small-/x estimate (44), i.e., Q;* « 
\/2/|/i|. In the limit p — )• 00 the value of Qf should 
converge toward the corresponding value in the finite-size 
behavior of the homogeneous system with open boundary 
conditions, see App. B. Thus 



lim Q;* = — , 

p— >oo Zj 



(54) 



where / is the filling factor given in Eq. (35). 

The asymptotic values can also be computed using 
the LDA of the total particle number, cf. Eq. (40). Since 
the level crossings occur at the boundary between the 
regions with N — k and N = k + 1, asymptotically we 
have k « Ni^a. = c{fi)l, thus 



or = i/c(m), 



(55) 



whose numerical values coincide with the estimates of Q* 
obtained above, see Table 1. This is actually a further 



evidence, beside the results of Sect. Ill, that the LDA of 
the particle number is asymptotically exact. 



2. Trap-size dependence of the gap 

For any < 1 and power of the confining potential, 
the results for the gap show the behavior 



A = /-iAa((/') + 0(/-'), 
Aa(</>) =ai(0), 



(56) 



with t{(i)) given by Eq. (49), and the coefficient a depends 
on p and /z. We determine the single parameter a by 
fitting the peak values to the form 



^peak 



= \al~'^ -\-cl~ 



(57) 



(the term is absent just for this quantity); the fit 
quality is usually excellent for I > 10. We report some 
results for the constant a in Table I. Note that the values 
for p = 2 appear to approach the small- estimate a = 
\pl with increasing /u. The value of 4> at the peak of the 
gap converges as (/'peak = l/2-|-d~^-|-0(/~^) with c <C 1, 
of the order of 10~^. 

We note that the modulated trap-size dependence of 
the gap within the superfluid region |^| < 1 appears to 
be largely universal, being substantially independent of 
/U and p, apart from a trivial normalization. 

This is also confirmed by the finite-size behavior of 
the homogeneous XX chain corresponding to the p — >■ 00 
limit of the confining potential, where the model becomes 
equivalent to a homogeneous chain of size L = 21 with 
open boundary conditions (more precisely, the p — )• 00 
limit corresponds to a chain with L = 2[l\ + 1 when the 
center of the trap coincides with the middle site of the 
chain, and L — 2[l\ when the center is in the middle 
between two sites). The results of App. B show that the 
modulated trap-size behavior of the gap found at finite 
values of p persists in the p = 00 limit. 



3. The particle density 

The particle density (no) at the origin shows a mod- 
ulated asymptotic behavior as well, but with a period 
twice the period of A; therefore it is useful to define 



(2k) 



(2k+2) _ ,(2fe) 



for / 



(2fe) 



< / < 







,(2fe+2) 
'0 



(58) 







thus < ^ < 2 (note that either ^ = 0or^ = (^-|-l 
for large I). The results of the numerical diagonalization 
show that (no) satisfies the asymptotic behavior 



(no) - / = /-M„(,/.) 



(59) 
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2.0 2.5 



FIG. 10: The functions h and g, cf. Eq. (60), vs. Y = 
pp/(p+i)^ at odd peaks of A at ^ = 0, for p = 2 (above) 
and p = 4 (below). Results for p = 6 are of the same quality. 
h is real, since kp — n/2. g at peaks is compatible with 0. h 
at even peaks is opposite in sign to h at odd peaks. The data 
at different trap sizes clearly approach TSS functions when 
increasing I. Here and in the following, in the legends I is 
rounded to the nearest integer for the sake of presentation. 




FIG. 11: The functions h (solid lines) and g (dashed lines), cf. 
Eq. (60), vs. Y = l~^^^x at odd peaks of A for p = 2 and /i = 
cos7r/5 (/ = 1/5). g at peaks is compatible with 0. h at even 
peaks is opposite in sign to h at odd peaks. The "noise" more 
evident at smaller I is due to the ambiguity of disentangling a 
"fast" oscillation e^'*^'' from a "slow" osciUation h{Y)/\h{Y)\ 
in Eq. (60). 




where / is the filling factor of the homogeneous model 
without trap, cf. Eq. (35). Some results for the constant 
b are reported in Table I. Note that the 0{l^^) term has 
zero average over the period. 

An analogous behavior of the particle density in the 
middle of trap is found in the limit p — >■ (X), as shown by 
the results of App. B for the finite-size behavior of the 
homogeneous XX chain with open boundary conditions. 

The particle density as a function of the distance x 
from the middle of the trap turns out to behave as 

K) « Pida(X) + ri Re{h{Y, ^)e2^^-^ + g(Y, 4>)} 

x = x/i, y = xrf/(f+i), (60) 

where kp — tt/ — arccos/x, g is real, and terms sup- 
pressed by higher powers of are neglected, g and 
h are discontinuous at (f> — 1. Note that two scaling 
variables X and Y appear in the above equation, distin- 
guishing the scaling behavior of the two terms. Note also 
that the leading term depending on X is the LDA of the 



FIG. 12: The functions h and g, cf. Eq. (61), vs. Y = r'^^^x at 
peaks of A for p = 2 and fi = 0; his real, since kp = n/2. The 
data at different trap sizes clearly approach TSS functions 
when increasing I. 

particle density, cf. Eq. (37). Some results are plotted in 
Figs. 10 and 11. We find that g{Y) = at peaks of A 
(<^ = l/2, 3/2). 

4- The particle density correlation 

Results for p = 2, 4, and 6 show also that the connected 
density correlation scales as 

G,,{x, 0) « r2p/(P+i) Re{h{Y, 4>)e^'''-^ + ~g{Y, 0)}, (61) 

where Y — a;Z~P/'-P+^\ and g is real, h and g are dis- 
continuous at ^ = 1. Results are shown in Figs. 12 and 
13 at = and p — 2, 4 respectively, and in Fig. 14 at 
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FIG. 13: The functions h (solid lines) and g (dashed lines), 
cf. Eq. (61), vs. Y = r^/'f+^'a; at = and at odd peaks of M = 
A for p = 4. Results for p = 6 are of the same quality. 



FIG. 15: x^^'^Gb{x,Q) at peaks of A vs. x/l for p = 2 and 
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FIG. 14: The functions h (solid lines: amplitude; dashed lines: 
phase/ (207r)) and g (dotted lines) vs. Y — l~^^^x at odd peaks 
of A for p = 2 and ^ — cos 7r/5 (/ = 1/5). 



fj, = cos7r/5 and p = 2, at the peaks of the gap. We find 
that g{Y) = \h{Y)\ at peaks of A (0 = 1/2, 3/2). 



5. The one-particle density matrix 

Another interesting quantity is the one-particle density 
matrix (30). In the homogeneous system without trap 

Gbixi^Xj) \xi - Xj\~^^'^, (62) 

from which we can read the RG dimension of the bosonic 
operator b in the superfluid phase, i.e., yb = 1/4. This 
behavior is also observed in the presence of the trap for 
sufficiently small distances [18]. Fig. 15 shows results for 
p = 2 and /x = 0: the data of x^^^Gb^x, 0) for different Is 
appear to collapse to a unique curve when plotted versus 
x/l, apart from small oscillations with a decreasing am- 
plitude for Z — > cx). Only the oscillations depend on 0. 



For X > I, Gb{x,0) is exponentially suppressed. These 
results indicate an asymptotic scaling behavior given by 

Gb{x,0)^r^^^g{x/l). (63) 

The region where Gb{x,0) appears to rapidly vanish, 
i.e., for x/l « 1 in Fig. 15, corresponds to the region 
where /XeS ~ 1, cf. Eq. (36), which is the value of the 
chemical potential corresponding to the superfluid to 
empty state transition, where the particle density of the 
ground state vanishes. We thus expect that, for generic 
values of /i andp, the region around x = Xc = 1(1 — iS)^^^ , 
where ^ceixc) — 1, develops critical modes related to a 
low-density Mott transition. The effective chemical po- 
tential can be expanded around 

Meff = fi+{x/ir = l + {l-^i)p^^^+0[{x-x,f]. (64) 

Thus, the behavior around Xc is essentially analogous to 
that arising at the low-density Mott transition n = 1 in 
the presence of a linear potential V/ ~ r/l. Around Xc, 
critical modes should appear with length scale ^ ^ 1°', 
where a is the exponent associated with a linear exter- 
nal potential. The value of a can be inferred by RG 
arguments analogous to those leading the determination 
of the trap exponent 6 at the low-density Mott transi- 
tion [23], which give a = 1/3.^ We thus expect that the 
transition region around x = Xc enlarges as 

Ax - (65) 



^ The exponent a can be determined by a RG analysis of the per- 
turbation corresponding to a linear potential Vi{x) = ux, i.e., 
J d'^ X dt Vi{x)\{j>{x)\^ , at the fixed point of the continuous theory 
describing the Mott transition [11]. The exponent cr is related to 
the RG dimension j/u of the parameter u, which can be obtained 
from the relations y-u ~ 1 = d + z — y\^\2 = j/^j = 2, thus yu = 3, 
and therefore cr = 1/yu = 1/3 for d = I and d = 2. 
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FIG. 16: Ax, defined by Eq. (66), vs. « for p = 2, 4 and 
fi = 0. The full lines shows a fit to Ax — al" leaving free 
the exponent a. The resulting estimates of a axe in perfect 
agreement with the RG prediction ct = 1/3 independently of 



independently of the power-law p of the confining poten- 
tial. We study this phenomenon numerically by comput- 
ing Gb(a;,0) for Is corresponding to odd peaks of A; we 
take 



Ax = l 



(66) 



where Xmax is the abscissa of the rightmost maximum of 
x^/'^Gh{x,Q)] the results agree very well with Ax ~ 1^^^, 
see Fig. 16. 



D. Quantum entanglement in the superfluid phase 

We divide the chain in two parts of length Ia and 
L — Ia, and consider the von Neumann entanglement en- 
tropy S{Ia] L) for open boundary conditions. In the ab- 
sence of the trap and for open boundary conditions, the 
von Neumann entanglement entropy is essentially deter- 
mined by the conformal field theory which describes its 
continuum limit [32], 



S{Ia;L)^ -\n[Lsin{nlA/L)] 
6 



(67) 



where c = 1 is the central charge corresponding to a 
relativistic free boson field theory. The 0{1) term, E{ii), 
depends on /.t; it is exactly known, see Eq. (Bll) (notably, 
E{id) = 0.287769699...). We consider specifically the half- 
lattice von Neumann entanglement entropy, 



S(L/2:L) = -\nL 
6 



E{f,) + Oil/L). 



(68) 



As shown in App. B, the amplitude of the 0{1/L) cor- 
rection is generally modulated by a function of ^ = 
2{[{L + l)f + l]/2} (where {x} = x - [x\ is the frac- 
tional part of x), giving rise to peculiar oscillations. 



In the presence of a trapping potential, and for L 
our numerical results show the behavior 



oo. 



Sl/2 



lim S{L/2-L) (69) 
\ hl[ae{^i) I] + E{fi) + Aoifi, + 0(1/P) 



for any p, where E{^) is the same constant in the absence 
of the trap, cf. Eq. (68). Eq. (69) defines an entanglement 
length scale [23] 



Ce = Oe I. 



(70) 



Some results for the amplitude Og are reported in Table 
1. For p — > oo we recover Eq. (68) with L = 21, thus 



lim ae{fJ.) 



(71) 



The amplitude Ao of the 0(1/1) term turns out to be a 
periodic function of the trap size, through the depen- 
dence on <j), i.e., the phase-like variable defined in 
Eq. (58). An analogous behavior is found in the limit 
p oo, i.e., in the case of a homogeneous system of size 
L with open boundary conditions, see App. B. 

We also mention that similar subleading oscillations 
are observed in the half-lattice entanglement entropy of 
the XX model with gradients [33], i.e., in the presence of 
a linear external field. 



E. Some notable relations 

We have already shown that in the gapless superfluid 
phase the asymptotic modulated power-law behavior of 
the gap and the particle density in the middle of the 
trap is largely universal, being independent of /i and p, 
apart from trivial normalizations. In the following we 
show that also the amplitudes for different values of p 
are strictly related; they can be derived from the p = oo 
limit which corresponds to the homogeneous system with 
open boundary conditions, see App. B. 

We note that, using the entanglement definition (70) 
of length scale, for any < 1 and any p including p — >■ 
oo, the asymptotic behavior of the gap and the particle 
density can be written as 



1 



(72) 
(73) 



with (j) and (j) defined in Eqs. (47) and (58) respectively. 
Indeed, one can check that, within the high accuracy of 
our numerical estimates, the results reported in Table I 
satisfy the relations 



aip;p) 
b{n;p) -■ 



1 



aeifJ-;p)' 



(74) 
(75) 
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Moreover, Eqs. (72) and (73) reproduce the results for 
the homogeneous system with open boundary conditions, 
see App. B, by replacing the entanglement length scale 

with L + 1, with and 4> given by the corresponding 
expressions (B2) and (B5) respectively. 

These results provide a strong numerical evidence of 
the following statement: in the superfluid phase the 
asymptotic trap-size dependence of smooth observables 
can be obtained by replacing L (or, more precisely L + l) 
with in the asymptotic behavior of the homogeneous 
system of size L with open boundary conditions. 

Another notable numerical relation is found at = 0: 
the numerical data of the half-lattice entanglement 51/2 
in the presence of the trap indicates that the entangle- 
ment length scale ^e, cf. Eq. (70), is exactly given by 



2Q^l, 



(76) 



for any p, where Qi is the asymptotic periodicity of the 
level crossings, cf. Eq. (53), which is exactly derived from 
the LDA of the total particle number, i.e., from the rela- 
tion 



dx arccosx^ 



(77) 



We are quite confident that Eq. (76) holds, since the 
numbers reported in Table I show that it is verified within 
the numerical accuracy of our estimates of the amplitude 
tte for p = 2 and p = 4, which is O(10~^). Moreover, it 
correctly reproduces the p — >■ oo limit = 21. 

In addition, we find that at /i = the behavior of the 
half-lattice entanglement S1/2 in the presence of the trap 
turns out to be consistent with the following formula 

^1/2 = I InCe + m + ^/o{4') + O(lAe) (78) 

for any p, where 

M,;^) = l i + for 0<^<1, 

\ -l-c(3/2-(/)) for 1<<^<2. ^ ^ 

with c w -0.150 for p = 2 and c « -0.044 for p = 4. 
Note that the behavior of the homogeneous system with 
open boundary condition, cf. Eq. (B15), is obtained by 
replacing ^ + 1 for the values cf) = 1/2, 3/2, which 
are the only possible values for the homogcineous system 
in a chain with even L, corresponding to odd and even 
L/2 respectively. 



F. Discussion 

The trap-size dependence of the half-lattice entangle- 
ment shows that the confining potential induces a length 
scale which behaves as ^ ^ Z for any power of the poten- 
tial, at least for smooth observables. In the framework 



of the TSS, this would imply that the trap exponent is 
9=1 independently of p. This value of 6 is also obtained 
from the trap-size dependence of the gap, which scales as 
A ~ apart from a periodic dependence on / of its am- 
plitude, cf. Eq. (56). Indeed, the exponent of the power 
behavior is expected to be z6, and in the superfluid re- 
gion we have 2; = 1. A consistent scaling is also observed 
in the case of the one-particle density matrix Gb {xi ,Xj), 
see Sec. IV C 5. On the other hand, the density correla- 
tions show clearly a coupling with modes at the Fermi 
momentum kp = tt/, which are apparently character- 
ized by a difi'erent length scale , scaling as with 

C = p/(p + i). 

These results may be explained by the nontrivial cou- 
pling of the confining potential with the free bosonic field 
4>{x) of the continuum theory, due to the fact that in 
the continuum limit the spin operator erf can be writ- 
ten as a sum of a slow contribution proportional to dxcj) 
and a rapidly oscillating contribution proportional to 
^2ikFX^^ see, e.g., Refs. [25, 26]. Assuming that the lat- 
ter term is suppressed for smooth or global quantities, 
such as the half-lattice entanglement entropy and the 
gap, where its effects should get averaged out, we may 
argue that 6* = 1 is indeed the expected trap exponent. 
A heuristic argument may be obtained by noting that 
the perturbation J dxV{x) {d(j){x)/dx) can be rewritten 
as J dx {dV (x) / dx) (p(x) by integration by parts, whose 
first-order perturbation vanishes because V{x) is even in 
x. We then expect that the leading contribution comes 
from next-to-leading terms, like J dxV{x) [d^ (t){x) / dx"^). 
RG arguments applied to this perturbation, taking into 
account that tj) is & free Bose field, lead to 9 = 1 in- 
dependently of p, which is the result emerging from the 
numerical diagonalization of the Hamiltonian. On the 
other hand, correlators which are nontrivially coupled to 
the modes at kp may show a different length scale due to 
the coupling of the confining potential with the staggered 
term in . This is indeed what we observe in the corre- 
lation of density operators which are directly related to 
the operator . 



V. MODULATED TSS AT THE n = 1 MOTT 
INSULATOR TO SUPERFLUID TRANSITION 

In the section we study the effects of the trap at the 
n = 1 Mott insulator to superfluid transition, i.e., at 
He = —1, where the filling factor of the homogeneous 
system is / = 1. The confining potential V{x) gives 
rise to a change of the particle density from (n) « 1 
in the middle of the trap to (n) = at large distance, 
passing through the gapless superfluid phase, see Fig. 9. 
Specifically, for h = —1, the particle density appears to 
vanish when ficS ^ 1, cf- Eq- (36), thus, x/l > 2^^"^. 

We recall that the behavior around = — 1 of the 
homogeneous BH model without trap is essentially anal- 
ogous to that at /X = 1, because of the invariance under 
the particle-hole exchange. At the n = 1 Mott insulator 



15 



I 

250 251 252 253 254 255 



0.4 




1 1 1 1 1 1 1 1 1 1 


- 


0.3 








0.2 






- 


0.1 








/ - 


0.0 


M 1 1 1* 




VI 1 1 IViVi 1 1 1/ 




0.3 




'I 1 1 1 1 


III 


1 1 1 1 




0.2 










- 


0.1 












0.0 


" , , , , 1 




1 1 1 1 


1 1 1 1 





250 251 252 253 254 255 

I 

FIG. 17: The rescaled energy gap i'^^A (above) and the 
rescaled particle density in the middle of the trap f {1— (no)) 
(below) vs. / for fi = —1 and p = 2, whose trap exponent is 
9 = 1/2. Results for p = 4 and p = 6 are similar. 
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FIG. 18: P^A (above) and l\l - (no)) (below) vs. 0, defined 
in Eq. (82), for p = 2 (S = 1/2). The extrapolation to I ^ oo 
is obtained by assuming 0{l^^^^) leading corrections. In the 
legends / is rounded to the nearest integer for the sake of 
presentation. 



to superfluid quantum transition, the critical exponents 
z and yfj_ and the trap-size exponent are the same as 
those at fi — 1, i.e., z = 2, = 2 and 9 = p/{2 +p). 
However, the particle-hole symmetry does not hold in 
the presence of the trapping potential, and the asymp- 
totic trap-size dependence appears more complicated at 
the n = \ Mott transition. This is essentially related 
to the presence of level crossings at finite values of the 
trap size, where the gap vanishes, as already found in 
the superfiuid region, for |^| < 1. As we shall see, the 
resulting trap-size dependence can be cast in the form of 
a modulated TSS, that is a TSS controlled by the same 
exponents as those at the low-density Mott transition, 
but modulated by periodic functions of the trap size. 



A. Modulated TSS of the gap 

Results for the gap and the particle density at the mid- 
dle of the trap are shown in Fig. 17. They suggest peri- 
odic asymptotic behaviors of the scaling quantities Z^^A 
and Z^(l — (no)) as functions of with a period given 
by the interval between two even (or odd) zeroes of the 
gap, and a marked difference between even- and odd- 
numbered crossings and peaks. 

In the large-Z limit we find that the interval between 



two even zeroes approaches a constant value, i.e.. 



P, 



(fc) 



(2/C+2) 



{2k) 



Pi* + 0(1' 



I — '■a '■0 ^ ^ I 

The large-Z limit is estimated to be 

P; ^ 1.11072073 for p = 2, 
Pi = 1.10243940 for p = 4, 
P; = 1.08184087 for p = 6. 



(80) 



(81) 



Analogously to the trap-size dependence in the superfluid 
region, see Sec. IV D, the asymptotic interval P;* can be 
estimated using the LDA of the total particle number, cf. 
Eq. (40). We find again identical results, showing that 
the LDA of the total particle density is asymptotically 
exact in the large-Z limit also at the n — 1 Mott tran- 
sition. The p — > oo limit of Pj* can be easily computed 
using the LDA, obtaining P;* — )• 1 for p — )■ oo. Correc- 
tions to the LDA of the particle density are discussed 
below. 

The asymptotic periodic properties clearly emerge 
from the results shown in Figs. 18, 19 and 20, where 
P^A and — (no)), for p = 2, 4, 6 respectively, are 
plotted versus the phase-like variable 



I -I 



{2k) 



,{2k+2) 



,(2fc) ' 



{2k) 



<l <l 



{2k+2) 
' 



(82) 
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FIG. 19: P^A (above) and ^^(1 ^ 
in Eq. (82), for p = 4 (S = 2/3). 



(tiq)) (below) vs. (j>, defined 




FIG. 20: (above) and f {I ^ 
in Eq. (82), for p = 6 (S = 3/4). 



(no)) (below) vs. (j!>, defined 



thus < (/)< 1. 

The results for the gap show that the quantity P^A 
approaches an asymptotic periodic function Aa(<?!>) in the 
large-^ limit. Therefore, they provide a clear evidence for 
a modulated asymptotic behavior 



^A(0)r^''[i + o(r«)] 



(83) 



with K w 9/2 (see below). Note that, for p — 2, the 
scaling of the gap is controlled by the same exponent 
found in the region < 1, because we have 20 = 1; on 
the other hand, for p — A and p = 6 we have respectively 
29 — 4/3 and 29 = 3/2, which are easily differentiated 
from 1. 

Wc should note that in these calculations the trap of 
size I is centered in the middle of the chain of size L; more 
precisely we consider odd Ls to have the center of the 
trap coincide with the middle site of the chain. Even Ls 
would instead correspond to traps centered between the 
two central sites. Unlike the previous cases, at /i = — 1 
this difference must be taken into account, but it gives 
only rise to an interchange of the role of the even and 
odd zeroes of the gap. Therefore, in the case of even Ls, 
one may simply redefine the phase-like variable as 



I -I 



(2fe-l) 




(2fc+l) ,(2fe-l) ' 



(2fe-l) 



<l <l 



(2fe+l) 




(84) 







The two different definitions (82) and (84) remove the de- 
pendence on the parity of L (cf. App. A) from the modu- 
lation functions A((j)) of the asymptotic behaviors of the 
observables, such as Aa(0)- Using, e.g., definition (82) 
for ah Ls would lead to A(0)|cvon l = ^({0 + 0o})|odd l, 
where {x} = x — [x\ is the fractional part of x. Note 
that centering the trap at an arbitrary distance q from 
the nearest site of the chain would modify significantly 
the modulation. 

By definition, ^a(O) = ^a(I) — 0. Aa has another 
zero 00 corresponding to the odd level crossings. Com- 
putations of (j)o for the XX model and p = 2, 4, 6 are 
shown in Fig. 21, where they are plotted vs L~^. They 
asymptotically approaches the value (j>o = 0.750000 for 
p = 2, 00 = 0.833333 for p = 4, 0o = 0.875000 for p = 6, 
and (po = 0.916667 for p — 10, with 0{l~^^) corrections. 
Note that these values are compatible with the simple 
formula 



0o = (p+l)/(p + 2). 



(85) 



Other features of the modulation function Aa(0), such 
as the location of the large and small peaks, (pi and 
4>s respectively, and its peak values Aa(0;,s), are ap- 
proached with power-law scaling corrections 0{l~'^). We 
find K K, 1/4 for p = 2, K « 1/3 for p — A, and k « 3/8 
for p = 6, which are consistent with k — 9/2. See, e.g.. 
Figs. 22 and 23. Note that k = 9/2 should be con- 
sidered as a phenomenological result, because we do not 
have theoretical arguments to derive it. 

A faster approach to scaling is found for the ratio 
A/ / As between subsequent large and small peaks of the 
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FIG. 21: (po vs. r*, for the XX model (6) at ^ = -1 and 
the XXZ model (33) at = -1/2 and = -1/2, for p = 2 
(above), p = 4 (middle), and p = 6 (below). 



gap, which is expected to be universal (essentially be- 
cause it is independent of normalizations), see Fig. 24. 
(We compute the ratio at Is corresponding to , obtain- 
ing A/ by cubic spline interpolation.) We find a behavior 
compatible with 
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FIG. 22: Scaling of the large (odd) and small (even) peaks of 
the gap and their location (pi and (ps for p — 2, for the XX 
model (6) and the XXZ model (33) at > = -1/2. 



It is worth noting that, in the case p — 2, A^{(j)) has 
apparently a simple form: it is approximately formed 
by two reflected similar triangles, as shown in Fig. 18. 
Indeed, this shape would require A//As = (j)Q/{l — cpo) = 
3, while we obtain a very close value Ai/As — 3.013 from 
the I — )■ oo extrapolation of the results at fixed I. 

The behavior in the large-p limit can be guessed from 
the results obtained at finite p shown above. Note that 
they do not apparently approach the finite-size behavior 
of a homogeneous BH model on a chain of size L — 21 
with open boundary conditions, as one may naively ex- 
pect. For example, the gap of the homogeneous BH 
model at /i = — 1 behaves as A = 0(1/L^) without show- 
ing level crossings. As a consequence, the phase-like 
variable does not have any corresponding quantity in the 
BH model without trap. Thus the relations between the 
large-p limit and the finite-size behavior of the homoge- 
neous BH model are not straightforward. 

Finally, let us mention that the gap at /i < — 1 shows 
nontrivial trap-size dependence as well, because the phe- 
nomenon of the level crossings persists. The numerical 



^ This result can be easily derived using the particle-hole exchange 
symmetry, and the corresponding result at the low-density Mott 
transition, at /i = 1, see, e.g., Eq. (19). 
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FIG. 23: Scaling of the large (odd) and small (even) peaks of 
the gap and their location (pi and (jts for p — 4, for the XX 
model (6) and the XXZ model (33) at = -1/2. 



results at fi = —1.2 show again a periodic structure of 
the gap, which gets suppressed as I /I for any p. They 
show large and small peaks, but the smallest one gets 
rapidly suppressed in the large-^ limit. We observe that 
the large peak turns out to scale as A; ~ 1/Hor p = 2 and 
p = 4, while the small peak decreases as ~ exp(— a^); 
therefore, the ratio A^/A; vanishes as I — > oo; in the 
same limit, ^ 1. Although for /i < — 1 the homo- 
geneous system without trap has a gap proportional to 
fis = — M ~ Ij see, e.g., Ref. [25], here we find that the 
periodic trap-size dependence of the gap tends to be sup- 
pressed as l/l in the large- Z limit. This is related to the 
fact that, in the presence of the trap, including its large-Z 
limit, we have always some critical regions, for example 
the super fluid vegiows between the {ux) ~ 1 and {rix) = 
plateaus, for 
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FIG. 24: Ai/As vs. for the XX model and for the XXZ 

model at jz = —1/2, for p = 2 (above), p = 4 (middle) and 
p = 6 (below). 



(-l-A^)'/"<N/^<(l-M)'^^ (87) 



which becomes larger and larger with increasing I. An 
analogous behavior is found for the trap-size dependence 
of the XY chain when the middle of the trap is in the 
quantum ferromagnetic phase [23] . 



B. The particle density and its correlators 

Concerning the particle density, we recall that (rix) = 1 
in the absence of the trap and the RG dimension of the 
particle density is y„ = 1 at the Mott transition. Thus, 
TSS predicts the scaling behavior (no) — 1 ~ for the 
particle density at the middle of the trap and = — 1. 
Like the gap, we find that the trap-size dependence of 



19 



0.25 



= 0.20 



0.15 



0.10 



- XX, even peak 
= XX, odd peak 
A XXZ, even peak 
V XXZ, odd peak 



0.0 



0.1 



0.2 



0.3 



-1/4 



0.9- 



Q^^-^r ° XX, even peak 

^ - = XX, odd peak 

- A XXZ, even peak 

7 - ' XXZ, odd peak 



0.6- 
0.0 



0.1 



V V V V 

' o!2 ' 



0.4 



0.5 



0.3 



I ^ 

0.4 



0.5 



D 



D 0.0 




FIG. 25: Scaling of the particle density at the peaks of the 
gap for p = 2 (above) and p = 4 (below), for the XX model 
and the XXZ model (33) at = -1/2. 



FIG. 26: The function 7) vs. X = I at odd and even peaks 
of A at /i = —1, for p = 2 (above) and p — 4 (below). 



(no) is described by a modulated TSS, i.e., 

(no) - 1 « r'Voi^b), (88) 

where </) is the same phase-hke variable defined for the 
scaling of the gap, of. Eq. (82). Figs. 18, 19 and 20 show 
1^(1 — (no)) vs. <j) ioT p = 2, 4, 6 respectively. The ap- 
proach to the asymptotic behavior is apparently charac- 
terized by power- law 0{l~'^) corrections with k ~ 6/2, 
analogously to the behavior of the gap. Notice that, un- 
like the superfluid case, the leading term has a nonzero 
average over the period. Results for the particle density 
at the peaks of the gap are shown in Figs. 25 for p = 2 
and p = 4. 

The spatial dependence of the particle density at large 
trap size turns out to be described by the following scal- 
ing behavior 

K) « pid^ix/l) + l-'ViX, cp), X = l-'^x. (89) 

As already found in the superfluid region, the term de- 
pending on x/l is given by the LDA of the particle den- 
sity, cf. Eq. (37). It plays the role of an analytical con- 
tribution which must be subtracted to observe scaling in 



the expectation value of the density operators at phase 
transitions [23, 24]. Some results for the scaling function 
I?(X, (p) at the peaks of the gap are shown in Fig. 26. 

As already shown by the data at /i = — 1 reported in 
Fig. 9, the particle density is quantitatively dominated 
by its LDA at large trap size, which scales as x/l. 

Assuming a modulated TSS, we expect that the den- 
sity correlation behaves as 

Gn{x) = {nmx)c - r^^GniX^dp). (90) 

This is confirmed by the numerical results. For exam- 
ple. Fig. 27 shows P^Gn{x) vs. X ior p — 2 and several 
values of I corresponding to peaks of A (i.e., (p = (pi or 
(ps asymptotically). The approach to scaling is clearly 
observed. 



C. The one-particle density matrix 

The modulated TSS also applies to the critical behav- 
ior of the one-particle density matrix. We recall that 
the RG dimension of the bosonic field bi is yi, — \/2 ai 
the Mott transitions. Consistently, the numerical results 
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FIG. 27: The rescaled connected correlation l^^G„{x) vs. X = 
xl~^ for p — 2, fi = —1, and Is corresponding to even and odd 
peaks of A. The data at different trap sizes clearly approach 
a TSS function when increasing I. 
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FIG. 28; The rescaled one-particle density matrix fGb{x) vs. 
X = xl^^ for p — 2 (above) and p = 4 (below), at /i = — 1, 
and for /s corresponding to odd and even peaks of A. 



FIG. 29: l^''^Gb{x,Q) vs. x/l for p 
taken at the odd peaks of the gap. 



show the behavior 



2 and n = —1. Data 



Gb(x, 0) EE {bib,) = {a^at) « r'MiX, <f). (91) 

In Figs. 28 we plot l^Gb{x, 0) vs. X for p ^ 2, 4, ^ = -1, 
and several values of I corresponding to peaks of A (i.e., 
<j} = <j)i 01 (j) = (j)s asymptotically). 

More generally, we may extend the modulated TSS 
including the dependence on /i = /i + 1 around fi — Q. 
Setting = P^' 



"/i, we expect 
G6(x,0) « r^7W(/i„X,, 



(92) 



Note that the modulated TSS is observed at fixed 
X = x/l^ , thus, since 9 < 1, the region oi x ~ I gets 
hidden in the limit X ^ 00. Fig. 29 shows numerical 
results for p = 2 versus x/l, where we clearly observe 
that the behavior for x ^ I \s approximately scaling as 
x/l. This is essentially related to the simple scaling of 
the effective chemical potential, cf. Eq. (36). However, 
a more careful analysis shows another interesting scaling 
behavior around the region 



= 2i/fZ, 



(93) 



where the spatial dependence of the effective chemical 
gives rise to a transition from the superfluid phase to the 
n — Q Mott phase. In this region the critical behavior 
should be governed by the linearized potential at Xc, cf. 
Eq. (64) , and therefore by the corresponding RG scaling. 
Analogously to what observed at = 0, see Subs. IV C 5, 
the transition region around Xc is expected to enlarge as 
Aa; ^ V-/^ independently of p. This RG prediction is 
fully supported by the same numerical analysis outlined 
at the end of Subs. IV C 5, for several values of p. 



D. The von Neumann entanglement entropy 



Finally, we consider the von Neumann entanglement 
entropy, and, specifically, the half-lattice entanglement 
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FIG. 30: The half-lattice von Neumann entanglement entropy 
Si/2 at /i = —1, for p = 2, 4, 6, 10. For p — 4, p — 6, and 
p = 10, the plotted data are already very close to the large-/ 
limit. 



entropy defined in Eq. (32) for chains with even L. We re- 
call that the half-lattice entanglement entropy S{L/2; L) 
of the homogeneous BH model vanishes at /i = —1, as 
it does at the low-density Mott transition point ^ — I, 
where it also vanishes in the presence of the trap, see 
Sec. II B 5. We instead find that at the n — 1 Mott tran- 
sition point fi = —1 the presence of the trap gives rise 
to nonzero values depending of the phase-like variable 0, 
i.e., 



S,/,^Asi4>) + 0{l-), 



(94) 



with K roughly compatible with 9/2. Results for p = 
2, 4, 6, 10 are shown in Fig. 30. Notice that the discon- 
tinuity in As occurs at = 0o = {p + ^)/{p + 2), cf. 
Eq. (85). The data suggest As {4>) ^ as p —> oo (for 
(/) < 1), consistently with the "naive" p — >■ oo limit of the 
homogeneous model with open boundary conditions. 



E. Universality of the modulated TSS 

We have shown that the BH model (1) is characterized 
by a modulated TSS at the n = 1 Mott transition. A 
natural question concerns its universality, i.e., whether, 
beside the critical exponents, also the modulation is uni- 
versal. 

In order to investigate this issue, we again consider 
the XXZ model (33), and present results obtained by 
DMRG calculations for = — 1/2, /i = /ic = —1/2 which 
corresponds to the n — 1 Mott transitions, and various 
power laws of the confining potential. The values of the 
trap size we could reach, I < 300, are significantly smaller 
than those for the XX chain, which are O(IO^). 

The numerical results show periodic structures of the 
gap and the particle density analogous to those of the 
XX chain, see Fig. 17. We find again level crossings 
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FIG. 31: The rescaled connected correlation l'^^Gn{x) (above) 
and the rescaled one-particle density matrix fGb{x) (below) 
vs. X = xl~^ for the XXZ model at jz ~ —1/2 and fi — fic, 
p — 4: and Is corresponding to peaks of A. For comparison, 
we also show results for the XX model for jj. = —1, p — 4 and 
a large trap size, see also Sec. VB and C. 



whose interval tends to a constant in the asymptotic 
large trap regime. The asymptotic interval between even 
level crossings, cf. Eq. (80), is given by Pj* — 1.550 for 
p = 2, P* = 1.293 for p = 4, and P^ = 1.200 for p = 6. 
Note that these values differ from those found for the XX 
model, see Eq. (81). 

After introducing a variable (f) defined as in Eq. (82), 
the gap turns out to behave analogously to the XX model, 
see Eq. (83). Results for the position (/)o of the odd zero 
are shown in Fig. 21: they are consistent with the asymp- 
totic values found for the XX model, although they show 
a slower approach, which is likely 0{l~^) (for the XX 
model it was 0{l~^^) for this particular quantity). Re- 
sults for the large and small peaks of the gap, and their 
locations, are shown in Figs. 22 and 23, respectively for 
p = 2 and p = 4. They appear to converge toward the 
same asymptotic values of the XX model, with similar 
corrections. The ratio Aj/Ag, which is inde- 
pendent of normalizations like (j)o, appears to scale ac- 
cording to Eq. (86), and the large-Z extrapolated values 
are consistent with those found for the XX chain, see Fig. 
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24. The agreement for p = 4 and p = 6 is very good; for 
p = 2, the data are clearly approaching the XX value, 
but the values of I considered are not sufficiently large to 
provide a precise extrapolation. 

The DMRG results for the particle density at the ori- 
gin, are consistent with the scaling behavior (88) with 
a universal scaling function ('/'), as also shown by the 
results at the peaks of the gap reported in Fig. 25. 

The behaviors of the density-density correlations func- 
tion and of the one-particle density matrix also support 
the universality of the modulated TSS. Examples are 
shown in Fig. 31 for p = 4. With increasing I the data 
appear to approach, although slowly, the same scaling 
curves found for the XX model. Note that the asymp- 
totic curves of different models are expected to be uni- 
versal apart from normalizations. This implies that they 
should match after an appropriate rescaling of the axes. 
Assuming universality, the results of Figs. 31 indicate 
that such rescalings are quite small for p = 4. 

In conclusion, the DMRG results for the XXZ model 
at jz = — 1/2 show a modulated TSS analogous to that 
observed for the XX model, thus supporting universality. 



VI. CONCLUSIONS 

In this paper we have studied the quantum critical be- 
haviors of confined particle systems described by the ID 
BH model (1) in the presence of a power- law confining 
potential V{r) = {r/iy, at the transitions between the 
Mott and superfluid phases, and within the gapless su- 
perfluid phase. We have considered the hard-core limit, 

— )■ oo, of the ID BH model, which allows us to study 
the effects of the confining potential by exact and very 
accurate numerical results. We have analyzed the vari- 
ous zero-temperature quantum critical behaviors within 
the TSS framework [23, 24]. In the following we list our 
main results. 

(i) At the low-density Mott transition, the TSS limit 
c;an be analytically derived within the quadratic spin- 
less fermion representation, leading to a Schrodinger-like 
equation for the lowest states. The RG scaling arguments 
leading to the TSS Ansatz (3) and (4), with critical ex- 
ponents z = 2, V = 1/2, and = p/{p + 2), are fully 
confirmed. We have computed the TSS of several observ- 
ables extending the results presented in Ref. [23], and 
checked universality by DMRG calculations within the 
XXZ model (33), which corresponds to adding nearest- 
neighbor density-density interactions in the BH model 
(1). The TSS functions show peculiar behaviors, whose 
main features, like the discontinuities in the scaling parti- 
cle density (see, e.g.. Fig. 1 for p = 2), are clearly related 
to the quantum nature of the transition. 

(ii) The trap-size dependence turns out to be more 
subtle in the region of parameters where the filling / of 
the corresponding homogeneous system is nonzero. This 
is essentially due to the presence of an infinite num- 
ber of level crossings of the lowest states when increas- 



ing the trap size. Nevertheless, the particle density of 
the ID hard-core BH model approaches its LDA in the 
large-/ limit, i.e., the value of the particle density of the 
homogeneous system at the effective chemical potential 
/Xeff(a;) = iJ,+ {x/l)P. Corrections are suppressed by pow- 
ers of the trap size, and show a nontrivial scaling behav- 
ior. 

(iii) The level-crossing phenomenon gives rise to a new 
interesting scenario at the n = 1 Mott transition, requir- 
ing a revision of the simplest TSS Ansatz (3) and (4) 
into a modulated TSS: the TSS is still controlled by the 
trap-size exponent 9 = p/{p + 2), as in the case of the 
low-density Mott transition, but it gets modulated by pe- 
riodic functions of the trap size. Indeed, the gap turns 
out to behave as A « Aa((^)Z"^^[1 + 0(/~^/^)], where the 
amplitude A/^{(j)) is a periodic function of the trap size I, 
through the phase-like variable 4> measuring the distance 
from the closest even level crossing, see Figs. 18, 19 and 
20 for p = 2, 4, 6 respectively. Modulated TSS is also 
observed for other observables, like the particle density 
and its correlators, and the one-particle density matrix. 
For example, the particle density {ux) at a distance x 
from the middle of the trap shows the asymptotic be- 
havior (rix) ~ p\da{x/l) + l~^'D{X,(j)) where X = l~^x, 
and p\da{x/l) is the LDA of the particle density, cf. Eq. 
(37). DMRG computations for the XXZ model (33) at 
jz = —1/2 and at the n = 1 Mott transition show an 
analogous modulated TSS, supporting its universality. 
The modulated TSS shows another peculiar aspect of the 
quantum nature of the Mott transitions. 

(iv) We have also studied the trap-size dependence in 
the gapless superfluid phase, whose corresponding con- 
tinuum theory is a conformal field theory with z = 1. 
In this region the asymptotic trap-size dependence turns 
out to be characterized by two length scales with dif- 
ferent power-law divergence in the large trap-size limit. 
One of them scales as ^ ^ Z and describes the behavior 
of observables related to smooth modes, such as the half- 
lattice entanglement; the other one scales as ^ ~ with 
C, = p/{p+l) and it is found in observables involving the 
modes at the Fermi scale kp = nf, where / is the filling 
of the homogeneous system. Moreover, the asymptotic 
power law behaviors appear modulated by periodic func- 
tions of the trap size. For example, the gap behaves as 
A ~ t{(j))l^^ where t((/)) is the triangle fmiction (49) and 
is a phase-like variable measuring the distance from the 
periodic level crossings, cf. Eq. (47). Some notable rela- 
tions are found for the behavior of observables related to 
smooth modes, such us the gap, the density at the origin, 
and the half-lattice von Neumann entanglement entropy. 
Their asymptotic behavior in the presence of a confining 
potential can be derived from that of the homogeneous 
BH model with open boundary conditions which repre- 
sents the p — > oo limit, by replacing the lattice size L 
with the entanglement length scale defined from the von 
Neumann entanglement entropy. 

The main features of the trap-size dependence reported 
in this paper should not be restricted to the hard-core 
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limit, i.e., the limit ?7 — > oo in the BH Hamiltonian (1). 
Specifically, the phenomenon of the level crossings at fi- 
nite trap size should persists at finite values of U, because 
the total particle number is conserved by the confining 
potential even at finite U. A study of the trap-size depen- 
dence at the Mott transitions of the ID BH model at fi- 
nite U would be important to further assess the extension 
of the universality of the modulated TSS observed in the 
hard-core limit. The presence of level crossings should 
also characterize the trap-size dependence in higher di- 
mensions. Therefore, modulated TSS may be also found 
at higher dimensions. 

Helpful discussions with P. Calabrese and M. Mintchev 
are gratefully acknowledged. 

Appendix A: Some details on numerical methods 

The numerical diagonalization of the Hamiltonian (8) 
is straightforward. Using lapack, L = 5000 chains can 
be diagonalized in a few minutes on a desktop PC. We 
consider chains of size L with open boundary conditions. 
The trap of size / is centered in the middle of the chain. 
The size L of the chain is taken sufficiently large to obtain 
results consistent with the infinite-size limit to machine 
precision. 

As it is clear from the discussion of the previous sec- 
tions, and specifically in Sec. HI, the density (ux) de- 
creases very rapidly for x > xida = 1 — fj)^^^: it is 
sufficient to choose L slightly larger than 2a;ida to have 
completely negligible finite-size effects. This allows us to 
obtain results for quite large O(IO^) trap sizes. A few 
examples of values of L giving boundary effects < 10~^^ 
for p = 2 are: jj, = 0.999, / = 10000: L = 1200; n = 0.9, 
I = 1000: L = 800; n = 0, I = 500: L = 1100; n = -1, 
I = 500: L = 1600. For larger values of p boundary ef- 
fects are smaller. Therefore, finite-L effects at fixed trap 
size are under complete control. [34] 

In the presence of the trap, most calculations are per- 
formed for odd Ls to have the trap centered in the middle 
site of the chain. We choose even Ls only to compute the 
trap-size dependence of the half-lattice von Neumann en- 
tanglement entropy. Notice that even and odd Ls may 
yield different results in the infinite-size limit, since the 
trap is centered between two chain sites and on a chain 
site respectively. For /i > — 1, this dependence on the 
parity on L disappears very rapidly with increasing I and 
it is totally negligible for / > 20 at /U = 0, for Z > 30 at 
H = -0.75 and for ? > 120 at /x = -0.9. For = -1, it 
must be taken into account, essentially because the role 
of the even and odd level crossings get interchanged, cf. 
Eqs. (82) and (84) in Sect. VA. 

Numerical results for the XXZ model are obtained by 
finite- volume DMRG; the Hamiltonian is not translation- 
invariant, therefore the initialization of the procedure 
is slightly nonstandard. The number of states M kept 
in the truncation procedure is chosen to have negligi- 
ble truncation errors; for the largest chains we consid- 



ered (L ^ 1000), we run with M up to 140, with dis- 
carded weights < 3xl0~^. Running at L = 999, p = 2, 
and /U = — 1 for a cycle of (j), determining 

^o^^\ /^ak ' and I'^^f^, cf. Subs. IV B, required about 100 
runs for different I, N combinations, each lasting about 8 
hours. 

Appendix B: Finite-size beliavior of tlie 
homogeneous ID hard-core BH model with open 
boundary conditions 

In this section we report some exact results for the 
finite-size behavior of the homogeneous ID hard-core BH 

model, or cquivalcntly for the homogeneous XX chain, 
with open boundary conditions. This formally corre- 
sponds to the p — )• oo limit of the confining potential, 
which becomes equivalent to a homogeneous chain of size 
L = 21 with open boundary conditions. More precisely, 
the p ^ 00 limit of the BH model with the trap corre- 
sponds to a chain with an odd L = 2[ZJ -I- 1 ([xj is the 
largest integer not greater than x) when the center of the 
trap coincides with the middle site of the chain, and to an 
even L = 2 [ZJ when the center is in the middle between 
two sites. 

In the infinite-size limit L — > oo, the filling / is given 
by [25] / = (rii) = (I/tt) arccos/x, thus in the range 
1 > fj, > —1 we have < / < 1. Let us now consider 
a homogeneous system of finite size L with open bound- 
ary conditions. The excitation number N for < 1 is 
exactly given by 

N=[{L + l)f\, (Bl) 

without any finite-L correction. For integer (L-|-l)/, the 
ground state is degenerate (A = 0); the lowest-energy si- 
multaneous eigenvectors of the Hamiltonian and the par- 
ticle number give N = {L + l)f and N = {L + l)f - 1. 

For / = 1/s with integer s, for every value of A'' we 
have a vacuum degeneracy when L -|- 1 = Ns, i.e., A = 
for L + 1 = (mods). For / = r/s with integer r and 
s, again A = 0forL-|-l = (mods), but we can satisfy 
L -|- 1 = Ns/r only for N = (modr). For irrational /, 
A never vanishes for integer L. 

Note that in the limit /i —1, thus / — > 1, Eq. (Bl) 
gives N — L without vacuum degeneration; this is the 
expected result for /x < —1. Analogously, Eq. (Bl) gives 
N — for ij, ^ 1 without vacuum degeneration, which is 
the expected result for ^ > 1. 

Eq. (Bl) suggests us to define 

(/)^{(L + 1)/}, (B2) 

where {x} = x — [x\ is the fractional part of x (i.e., the 
sawtooth function). For integer (L + 1)/, it is useful to 
label the two degenerate vacua with = 0,1 according 
toN + ^ = {L + l)f. 

For each value of /x, we observe that LA vs. (f) collapses 
on a curve proportional to the triangle function t(<p) de- 
fined in Eq. (49), with 0(1/L) corrections. Note that in 
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FIG. 32: The rescaled gap {1 - ^J.^)-'^/^^l-^{L + 1)A vs. <j), 
defined in Eq. (B2), for tlie liomogeneous system witii open 
boundary conditions, for several values of jj,. The apparent 
solid line is due to the accumulation of data points for fi = 
3/4, corresponding to an irrational /. Data for rational / are 
taken at L « 1000. 



the case of rational filling ip takes only a discrete set of 
values; for ^ = 0, e.g., takes the values and 1/2, cor- 
responding to odd and even L respectively. In Fig. 32 we 
plot data of {L + 1) A for several values of ^, correspond- 
ing to rational and irrational fillings (we use {L -\- 1)A 
rather than LA, which gives smaller, but comparable, 
finite-L corrections). Thus we have the asymptotic be- 
havior 



LA^aAt{(f>)+0{l/L), 



(B3) 



This corresponds to an asymptotic periodicity of the L- 
dependence of LA with period 1//. The numerical re- 
sults for OA turn out to be perfectly reproduced by the 
simple formula 



OA = 7r(l 



2n1/2 



(B4) 



Note that oa = at = ±1, indeed A = 0(1/L^) 
at /z = ±1 without level crossings, consistently with the 
fact that the corresponding continuum theory has z = 2. 
For 1^1 > 1 we instead have A = 2(|/^| - 1) + Oi^L-"^). 

We have also studied the particle density in the middle 
of the chain (rio), which is only defined for odd L. We 
define 



= 2{[(i + l)/ + l]/2}, 



< (/) < 2; note that either 
find that, at fixed /i. 



(i + l)((no)-/)-l-0. 



(B5) 
1 . Then we 

(B6) 



without any finite-i correction. We plot in Fig. 33 the 
r.h.s. of Eq. (B6) vs. for several values of /x; note that, 
for f = r/ s with integer r and s, is limited to a discrete 
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FIG. 33: The l.h.s. of Eq. (B6) vs. (f), defined in Eq. (B5), for 
odd Ls in the range 3 < L < 999. 



set of s values: {(2i — l)/s, i — 1, . . . , s} for odd s and 
{(2j — 2)/s, i = 1, . . . , s} for even s_(in this case, we also 
use = 2, which is equivalent to = 0, to distinguish 
the two vacua: = 2 for the N — {L + l)f vacuum and 
= for the N = {L + l)f — 1 vacuum. On the other 
hand, for irrational /s the values of are distributed all 
over the interval (0,2). 

For rational /, Eq. (B6) implies a very simple rational 
form for (np). For /i = 0, we have 



(L + l)((no) 



±1, L=l (mod 4), 
0, Lee 3 (mod 4), 



(B7) 



where the ± signs apply to the two degenerate lowest 
states with particle number N = {L + l)/2 and N = 
(L — l)/2 respectively. For ^ — 1/2, we have 



(L + l)((no)-i)=<^ \ 



0, 



L 
L 
L 



(mod 3), 

1 (mod 3), 

2 (mod 3). 



(B8) 



Eqs. (B7) and (B8) are indeed verified for all Ls to ma- 
chine precision. For other /is giving rational /, we obtain 
comparable results. 

For /i = 0, (rii) is given by a very simple expression for 
all points i ^ x + ^{L + 1), i = 1, . . . , L: 



(L + l)((n.) 



0, even L, 

0, odd L, even i, 

±1, odd L, odd i. 



(B9) 



where again the ± signs apply to the two degenerate 
lowest analogously to Eq. (B7). Eci. (B9) for « = ^(L + l) 
gives back Eq. (B7). 

Note that all the above formulae are invariant under 
the particle- hole exchange n,; — > 1 — n^, which implies 
N ^ L - N, f ^ 1- f, and fi^ -fi. 
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of the half-lattice behaves as 

S{L/2- L) = i ln(L + 1) + E{^l) + ^Ao{^l, 0) -f 0(l/i'), 

(B12) 

where <f) is defined in Eq. (B5), and j4o(/i, 0) has zero 
average over values of 0. Aq can be written as 



6 -|- ucf) 
-b-u{2- 



for 
for 



< 

1 < 



<1, 
< 2, 



(B13) 



where b and u depend only on /i. Indeed, in the case of 
irrational filling /, we find that, for sufhcicntly large L, 
the quantity 



FIG. 34: The function Ao{fJ.,<j)) as derived from data of the 
quantity reported in Eq. (B14), vs. (j), cf. Eq. (B5), for several 
values of fi. For = 3/4, all chains with even 200 < L < 2000 
are plotted. These results show that in all cases Ao can be 
parametrized as in Eq. (B13). 



We now consider the half-lattice von Neumann entan- 
glement entropy S{L/2;L) for even L and open bound- 
ary conditions in the supcrfluid phase, whose asymptotic 
large-L behavior is [32] 



1 



S{L/2-L) 



\nL + E{y) + 0{l/L). 



(BIO) 



The 0(1) term E{^) depends on /i; using the results 
reported in Refs. [28, 35], we derive 



dt 



iln(4/7r) 



cosh(t/2) 
4sinh^(V2) 



(Bll) 



2t sinh^(t/2) 



We recall that the von Neumann entanglement entropy 
S{L/2; L) vanishes at ^ = ±1. 

An accurate numerical analysis using data up to i = 
O(IO^) shows that, for any [yul < 1 including those corre- 
sponding to irrational filling, the entanglement entropy 



L[^(L/2;L)--ln(L 



1) 



Eiii)] (BI4) 
collapses on a single curve Ao(/i, 0) at fixed /i, given by 
Eq. (B13), with 0{1/L) corrections. In the case of ratio- 
nal /, we have a discrete set of possible values of 0, and 
the data accumulate at points located along two lines 
as described by Eq. (BI3). Some results are shown in 
Fig. 34. We find b ^ 1.48, u ^ -0.59 for /i = 3/4; 
b = 1.37, u ^ -0.52 for /i = 1/^2; 6 ^ 0.91, u ^ -0.17 
for /i = cos(27r/5) ^ 0.309017. 

We mention that subleading oscillations in the behav- 
ior of entanglement entropies have been also reported in 
other studies, see, e.g., Refs. [36, 37]. 

In the case ^ = 0, the values of corresponding to 
even L are restricted to (j) = 1/2, 3/2, which implies 
Ao ~ (—1)^/^. More precisely, our numerical results are 
accurately reproduced by the formula 

5(i/2;i) = iln(i + l) + i?(0)- (-1)^/2^^^^ 

+ 0{1/L'), (B15) 

with £:(0) = 0.28776969994598..., cf. Eq. (Bll). This ex- 
pression for the half-lattice entanglement entropy is con- 
sistent with an analogous formula reported in Ref. [36] for 
the XX model at ji = 0. The 0{1/L^) term remaining in 
Eq. (B15) is very small ('--' IQ^^L^^) and without oscilla- 
tions within numerical precision [^ 10~^^ on S{L/2; L)]. 
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